Quantum random walks - an introductory overview 
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This article aims to provide an introductory survey on quantum random walks. Starting from a 
physical effect to illustrate the main ideas we will introduce quantum random walks, review some 
of their properties and outline their striking differences to classical walks. We will touch upon 
both physical effects and computer science applications, introducing some of the main concepts 
and language of present day quantum information science in this context. We will mention recent 
developments in this new area and outline some open questions. 



God does not play dice. 
Albert Einstein 



I. OVERVIEW 

Ever since the discovery of quantum mechanics people 
have been puzzled by the counter-intuitive character of 
the laws of nature. Over time we have learned to accept 
more and more effects that are unimaginable in a classical 
Newtonian world. Modern technology exploits quantum 
effects both to our benefit and detriment - among the 
memorable examples we should cite laser technology and 
not omit the atomic bomb. 

In recent years interest in quantum information the- 
ory has been generated by the prospect of employing its 
laws to design devices of surprisingpower . New ideas 
include quantum cryptography 13, y] and quantum com- 
putation. In 1994 Peter Shor jj] discovered a quantum 
algorithm to factor numbers efficiently (that is in time 
that grows only polynomial with the length of the num- 
ber to be factored), which has unleashed a wave of activ- 
ity across a broad range of disciplines: physics, computer 
science, mathematics and engineering. 

This fruitful axis of research has uncovered many new 
effects that are strikingly different from their classical 
counterparts, both from the physical point of view as well 
as from a computer science and communication theory 
perspective. Over time these communities have gained a 
greater understanding of the concepts and notions of the 
other. The idea that information cannot be separated 
from the physical devise that is carrying it ('Information 
is physical') has firmly settled among them and has led 
to fascinating new insights. Acquaintance with the basic 
notions of each of these fields seems instrumental in the 
understanding of modern quantum information process- 
ing. 

In this article we will follow the trajectory of one of 
the many surprising aspects of quantum information; it is 
dedicated to quantum random walks. We will give a thor- 
ough introduction to the necessary terminology without 
overburdening the reader with unnecessary mathemat- 
ics. Starting with a very intuitive and physical example 
we will step by step introduce the language and notation 
of present day quantum information science. We will 



present the necessary background from computer science 
needed for a physicist to understand and appreciate the 
developments and results, assuming some rudimentary 
background of quantum mechanics, but no knowledge of 
computer science or quantum information theory. An 
excellent comprehensive introduction to quantum infor- 
mation and computation can be found in [l|. 

In this journey at the interface of several traditional 
disciplines we will follow quantum random walks both 
through physics and computer science. This tour through 
quantum information science will allow us to survey some 
of its most pertinent concepts and ideas today, like quan- 
tum algorithms, quantum computing machines, speed- 
ups, physical implementation, quantum circuits and de- 
coherence. It will be our mission to show how physi- 
cal phenomena translate into new computer science algo- 
rithms and vice versa. The idea is that a reader familiar 
with the standard idiom of quantum mechanics but un- 
accustomed to the language of 'qubits' and 'gates' will 
learn about quantum information science while learning 
about quantum random walks and their fascinating be- 
havior. We will survey some recent developments and 
results omitting most proofs but trying to develop an in- 
tuition. This article aims to bring the interested reader to 
a point where he or she can read and understand current 
research articles on the topic and develop an understand- 
ing for the interesting problems and open questions in the 
area. 

The structure of this article is as follows: First we give 
some physical intuition for random walks, providing a 
general flavor of the phenomenon in a physical setting 
(Section [nj. This is followed by a more rigorous defi- 
nition together with some necessary terminology to in- 
troduce the two main models of quantum random walks 
in Section IIIII We then present some computer science 
and probability background and mention some of the 
important algorithmic results stemming from quantum 
random walks in Section llVl Section switches back to 
physics and studies how these random walks could be im- 
plemented in some real physical system. Finally we deal 
with the more philosophical question of how the classical 
world emerges from quantum behavior via decohcrence 
using random walks as our example (Section IVI|) . All 
along the way we will outline open questions and future 
directions. 
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II. A GENTLE INTRODUCTION 

To illustrate quantum random walks and give an in- 
tuition for what follows we start with an example. This 
example - which might be thought of as a precursor to 
later models - is taken from the work of three physicists 
in 1993, Y. Aharonov, L. Davidovich and N. Zagury p|. 
Their work for the first time coins the term 'quantum 
random walk'. 

Imagine a particle on a line whose position is described 
by a wave-packet \ip Xo ) localized around a position xq, 
i.e. the function (x\ip Xo ) corresponds to a wave-packet 
centered around xq. Let P be the momentum operator. 
The translation of the particle, corresponding to a step 
of length I can be represented by the unitary operator 
Ui = exp(-iPl/K), so that Ui\tp xo ) = \ip xo -i). Now let 
us also assume that the particle has a spin-i degree of 
freedom. Let S z represent the operator corresponding to 
the z component of the spin and denote the eigenstates 
of S z by | t) and | |), so that S,| T) = || T) and S z \ | 
) = — 1|| I). From now on we will set K = 1 to simplify 
the notation. 

A spin-i particle is usually described by a 2- vector 

I*) = (IV^)i l^)) T > where the first part is the component 
of the wave-function of the particle in the spin- | T) space 
and the second one is the component in the | !)-space. 
Normalization requires that |||V^)|| 2 + |||V^)I| 2 — 1- To 
emphasize the tensor-structure of the space of the parti- 
cle we will write this in a slightly different but equivalent 
way as \^f) = aJ| |) ® \ipt) + a J -| J.}® where we nor- 
malize the two wave-functions (x/j^\tp') = (ip*-] ip*) = 1, so 
that |cJ| 2 + I ct-L 1 2 = 1. The tensor product '<g>' separates 
the two degrees of freedom, spin and space, and will allow 
us to view the resulting correlations between these two 
degrees of freedom more clearly. The time development 
corresponding to a translation by I on the larger state- 
space of the spin-i particle can now be described by the 
unitary operator U — exp(—2iS z <g> PI). This operator 
induces a kind of conditional translation of the particle 
depending on its internal spin-degree of freedom. In par- 
ticular if the spin of the particle is initially in the state 
I T)j so that its wave-function is of the form | |) (g) ^ 

then application of U transforms it to | t) ® I^I -j) and 
the particle will be shifted to the right by I. If the spin 
of the particle is in the state | |), i.e. the total wave- 
function is given by | !) <£> iV'i )' then the translation 
operator will transform it to | J.) ® IV'J +j) and the parti- 
cle will be shifted to the left. More interesting behavior 
occurs when the initial spin state of the particle, local- 
ized in xq, is not in an eigenstate of S z , but rather in a 
superposition 

|* i „) = (aT|T) + a i ||))®|^o)- (1) 

Application of the translation operator U will induce a 
superposition of positions 

U\V in ) = cJ| T) ® |^xo-i) WU) ® |V* +z>- (2) 



If at this point we decide to measure the spin in the S z ba- 
sis, the particle will be either in the state | |) ® \ip Xo _i), 
localized around xo + I with probability p^ = |cJ| 2 or 
in the state | J.) ® \tj) xo -i), localized around xq — I with 
probability p- 1 = |a^| 2 . This procedure corresponds to a 
(biased) random walk of a particle on the line: we can 
imagine that a biased coin with probabilities of head/tail 
p^ /p- 1 is flipped. Upon head the particle moves right and 
upon tail the particle moves left. After this step the par- 
ticle is on average displaced by I (jv —p^ ) . If we repeat this 
procedure T times (each time measuring the spin in the 
basis {| t)j I I)} and re-initializing the spin in the state 
«T| -|- Q.i| the particle will be displaced on aver- 
age by an amount (x) = T/(p T - p l ) = T/(|a T | 2 - la^ 2 ) 
and the variance of its distribution on the line will be 
cr 2 = 4in|cJ| 2 |cr>-| 2 = 4Tl 2 ptp±. This is exactly what we 
obtain if the particle performs a (biased) random walk 
on the line. We can imagine that the spin of the par- 
ticle takes the role of a coin. This coin is 'flipped' and 
the outcome determines the direction of the step in the 
random walk of the particle. 

Now let us consider the following modification of this 
procedure. Instead of measuring the spin in the eigenba- 
sis of S z we will measure it in some rotated basis, given 
by two orthogonal vectors {|s+), \s~)}. Alternatively we 
can rotate the spin by some angle 9 before measuring it 
in the S z eigenbasis. Let us be more formal and set up 
some more of the language used in quantum information 
theory. If we identify 

we can write 

Sz = l(o -1) = |d tXT I - [ i>U I)- (4) 

A rotation of the spin can be described by the matrix [53j 
m= (co S 9 -sin0\ 

Let us denote the measurement in the S z basis by M z . 
To explore the effect of the suite of operations M Z R(9)U 
on the initial state |^j n ) we have to slightly rewrite the 
operator U : 

jj = £ -2iS z ®Pl = g-idTXH-UXil)®^ 

= (\m\®e- iPl )(\l)il\®e iPl )- (6) 

The second equality follows from the fact that | |) (1 \®Pl 
and I |) (| I ® PI commute |54j together with the ob- 
servation that I t)(T I and | J.) (J, | are projectors, so 
I TXT | fe = I TXT I and I |)U | fc = I I) (i I when we 
expand the exponential. Now 

U\V in ) = (cJ| T ) e -Pi + a i I I) g, e lPl )\^ Xo ) (7) 

and applying R{9) from the left gives 

[(a T cos 9e~ lPl - sin 9e lPl ) \ |) ,„s 
+ (aUm9e- lPl + a 1 cos9e lPl )\ |)] ® \tp Xo ). [ ' 
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If the width Ax of the initial wave packet is much 
larger than the step length I we can approximate 



exp(±iPl)\tp x 



(/ ± iP/)l^xo 



With this in mind 
we can establish the state of the particle after the mea- 
surement M z : 



M z R(9)U\* in ) = 



with probabilities 



t)8(J-iPI^)|^) 
|) <g> (I - iPlS^)\ip xa ) 



(9) 



P 



T _ 



or cos 9 — or sin6>| 2 = \cv sin 9 + cos #1^10) 



and displacements 



W T := I 



a 1 cos 9 + a 1 sin ( 



IS 1 



a- 1 cost 



a cosf 



a* sin ( 



sin 9 + a-l cos 



If the width Ace of the initial wave packet is much 
larger than 16^ and 16^, we can again approximate 
{I-iPW^ X0 ) «exp(-iPWT.Jt)|^ o ) = |^ _ WT ,i). 

Now it is crucial to note that the displacement of 
the particle in one of the two cases, say, can be 
made much larger than I. For instance we may choose 

tan0 = IcJ/o^K 1 + e ) witn < N < !• Then the 

displacement of the particle in case we measured | f ) will 
be i^l i=a —2l/e which can be several orders of magnitude 
larger than I. However this event is 'rare': the proba- 
bility to measure | f) is p' ~ |a^a^| 2 e 2 . For the same 
choice of parameters the probability to measure | J.) is 



IS 1 



1 — \a' a 
i l(\a^ 2 - 



e and the displacement in that event is 
h Oile). For the average displace- 
ment p^W +p^l6* we get as before Z ( | ct^ | 2 _ |q^| 2 ), and 
similarly the variance of the walk will be unchanged. We 
have the effect that even though our translation opera- 
tion displaces only by I, in some (rare) cases the particle 
will jump much further than I. The effect is illustrated 
in Fig. n below. 

This is strikingly different from any classical behav- 
ior. Quantum mechanics allows us to post-select events 
that cannot be observed classically. Even though these 
events are rare events we cannot create a classical set-up 
that allows us to observe such an effect. This 'Gedanken- 
experimcnt' unveils some of the 'weirdness' of quantum 
mechanics. In what follows we will expand and show 
how variations on this phenomenon might be useful in 
the context of modern quantum information processing. 



III. THE QUANTUM RANDOM WALK 

Since the work by Aharonov et al. [5| mentioned in the 
previous section two models of quantum random walks 
have been suggested. They can be viewed as a gener- 
alization of the effect we have described above in one 
way or other, yet they differ in several aspects both from 
the example in Section [H] as well as from each other. 
In this section we will define both models and describe 
some results on the difference in behavior of classical and 
quantum walk. 



a) 
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FIG. 1: a) The particle is localized in position xo and enters 
the conditional displacement box U. Here depending on its 
spin degree of freedom it is shifted up or down by I. Upon exit 
its spin is measured in the z basis (Mz)- With probability 
1/2 it is found displaced by I resp. —I. b) The same set-up 
except that a rotation is applied to the spin after the particle 
exits the !7-box. After the measurement M z of the spin the 
particle now can be shifted up by far more than I (with small 
probability). With large probability it is shifted by less than 
I in the other direction. 



The key idea behind quantum random walks is to iter- 
ate the walk we have introduced in the previous section 
without resorting to intermediate measurements. In other 
words we will repeat the succession of unitary translation 
U and rotation R(9) without the measurement M z at in- 
termediate time steps. Furthermore we will discretize the 
position space of the particle, or in other words we put 
the particle on a lattice or a graph. Usually this space 
will be finite. A discrete and finite state space is impor- 
tant in the context of simulability on a (finite) computer. 
A quantum computer works with discrete registers, i.e. 
its state-space is a Hilbert space of large but finite di- 
mension. Discrctizing the quantum walk will allow us to 
map it to a computation by such a machine. 

Apart from the interesting new physics that we will 
uncover in the study of quantum random walks we will 
want to employ some of the 'weird' effects of the walk 
provided by quantum mechanics to enhance the comput- 
ing power at our disposal. The idea is that a quantum 
computer can implement (simulate) the quantum random 
walk efficiently and use this simulation to solve certain 
computational tasks. In this vein we hope to use prop- 
erties of quantum random walks to find more efficient 
algorithms on a quantum computer. 



A. The discrete model 

Let us define the formal set-up for a first model, the so 
called discrete time quantum random walk. Interestingly 
this model appears already in the works of Feynman [g 
with the discretization of the Dirac equation in mind. 



In the era of quantum information it was rediscovered 
in works by Meyer 0, 13 m connection with quantum 
cellular automata and by Watrous in the context of 
halting of quantum Turing machines and in a slightly 
modified version (with measurements) in [Tol ] thinking 
about space-bounded quantum computation. As a pos- 
sible computational tool and in a formal way it has been 
introduced and analysed in 2001 by Aharonov, Ambai- 
nis, Kempe and Vazirani (on general graphs) and by 
Ambainis et al. 01 ( on the line). 

We will define our model in one dimension, on the line 
or the circle. However the definitions carry over to the 
general case with slight modifications, which we will men- 
tion later. 

Let Hp be the Hilbert space spanned by the positions 
of the particle. For a line with grid-length 1 this space 
is spanned by basis states {\i) : i £ Z}; if we work on 
a circle of size N we have Hp = : i = . . . N — 1}. 
States in this basis will take the role of the wave function 
of Section |H] with \i) corresponding to a particle lo- 
calized in position i. We will not be concerned any more 
about the width of the distribution of the wave function 
- our model describes unitary transformations of states 
of a finite Hilbert space H and the notion of a particle 
is only used to guide our intuition from now on. Yet 
later, when we talk about the physical implementation 
of the quantum random walk, we will see that we can 
approximate the states of this mathematical model by 
particles (atoms, photons, etc.) with wave functions of 
finite width. 

The position Hilbert space Hp is augmented by a 
'coin'-space He spanned by two basis states {| |) , | J.)}, 
which take the role of the spin-^ space in the previ- 
ous section. States of the total system are in the space 
H = He <8 Hp as before. The conditional translation of 
the system (taking the role of exp(— 2iS z <g> PI)) can be 
described by the following unitary operation 

S = I TXT I ® E l i + + 1 Mi I ® E I*'" ^K 12 ) 

i i 

where the index i runs over Z in the case of a line or 
< i < N — lin the case of a circle. In the latter case we 
always identify and N, so all arithmetic is performed 
modulo N. S transforms the basis state | T) ® \i) to 
| t> ® |* + 1) and | |) ® \i) to | |) <g> |t - 1). 

The first step of the random walk is a rotation in the 
coin-space, which in analogy to the classical random walk 
we will call 'coin-flip' C. The unitary transformation C 
is very arbitrary and we can define a rich family of walks 
with different behavior by modifying C. For the moment 
we will want C to be 'unbiased' in the following sense. 
Assume we would initialize the walk in a localized state 
|0) with the coin in one of the basis states, say | T). If 
we chose to measure the coin-register of the walk in the 
(standard) basis {| |), | J,)} after one iteration (coin-flip C 
followed by translation S) we want to obtain the classical 
probability distribution of the unbiased walk, namely a 
translation to the right (| 1) ) with probability | and a step 
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FIG. 2: Galton's Board (Quincunx): 

The Quincunx is a device which allows a bead to drop through 
an array of pins stuck in a board. The pins are equally spaced 
in a number of rows and when the bead hits a pin it is equally 
likely to fall to the left or the right. It then lands on a pin 
in the next row where the process is repeated. After passing 
through all rows it is collected in a slot at the bottom. The 
distribution of beads approaches a Gaussian after many rows 
of pins. 

to the left (| — 1)) otherwise. A frequently used balanced 
unitary coin is the so called Hadamard coin H 

It is easy to see that the Hadamard coin is balanced: 

IT)®|0) -L(|0> + |1»®|0> 

^(|t)®|l> + |i>®|-l». (14) 

Measuring the coin state in the standard basis gives each 
of {| T) ® I I) ® | - 1)} with probability 1/2. Af- 
ter this measurement there is no correlation between the 
positions left. If we continued the quantum walk with 
such a measurement at each iteration we obtain the plain 
classical random walk on the line (or on the circle) . The 
distribution of this walk can be well illustrated with the 
help of Galton's board (also called Quincunx) in Fig. 
Its limiting distribution on the line (for large number 
of iterations T) approaches a Gaussian distribution with 
mean zero and variance a 1 = T. On the circle this dis- 
tribution approaches the uniform distribution (more on 
convergence of random walks on finite graphs in Sec. IIVII . 

In the quantum random walk we will of course not 
measure the coin register during intermediate iterations, 
but rather keep the quantum correlations between differ- 
ent positions and let them interfere in subsequent steps. 
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Asymptotic Pobability 



FIG. 3: The probability of being at position i after T steps 
of the classical random walk on the line starting in 0. 
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FIG. 4: The probability of being found at position i after 
T steps of the quantum random walk on the line, with the 
initial state = | J.) ® |0). Note that this distribution 

starts to differ from the classical distribution from T — 3 on. 
Furthermore the quantum random walk is asymmetric with a 
drift to the left. 



This interference will cause a radically different behav- 
ior of the quantum walk 55] . In particular we will see 
that the limiting distribution of the walk on the line will 
not approach a Gaussian and the variance a 2 will not be 
linear in the number of steps T. 

The quantum random walk of T steps is defined as the 
transformation U T , where U, acting on 7i = He <B> Hp is 
given by 



U = S-(C(g>I). 



(15) 



To illustrate the departure of the quantum walk from 
its classical counterpart let us evolve the walk (without 
intermediate measurements) for some steps starting in 
the initial state |<£>m) = | j) <£> |0) and study the induced 
probability distribution on the positions. 



|4>; 



u J_ 

1 

2 

u 1 



(|T>®|1>-Ii>®|-1» 

(lt>®|2)-(|t>-U»<8)|0) + |i)<g)|-2)) 

(I T)®|3) + |i)®|i) + |t)®|-i> 

(16) 



2V2 

-2\ I)® 1-1) -I i>® I -3)) 



This example already shows that the probability distri- 
bution induced by the quantum walk differs from the 
classical one. In Table|3]we list the classical distribution 
and in Table 0] the corresponding quantum distribution 
if we measure the position register after T steps. 




FIG. 5: The probability distribution of the quantum random 
walk with Hadamard coin starting in | [) CS> |0) after T — 100 
steps. Only the probability at the even points is plotted, since 
the odd points have probability zero. The dotted line gives a 
long-wavelength approximation (labeled P s i ow since it keeps 
only the slowly varying frequencies [15|), clearly showing the 
bi-modal character of the distribution. 



Figure [3] shows the probability distribution after T 
I i 10 steps of the quantum walk starting in | J.) ® |0). It is 
evident that the interference pattern of this walk is much 
more intricate than the Gaussian obtained in the classical 
case. One can clearly discern a two-peaked distribution. 

A first thing to notice is that the quantum random 
walk induces an asymmetric probability distribution on 
the positions, it is 'drifting' to the left. This asymmetry 
arises from the fact that the Hadamard coin treats the 
two directions | f) and | J.) differently; it multiplies the 
phase by —1 only in the case of | J,). Intuitively this in- 
duces more cancellations for paths going right-wards (de- 
structive interference), whereas particles moving to the 
left interfere constructively. There are two ways to mend 
this asymmetry. Inspecting Eq. (|14f) which describes the 
first step of the Hadamard walk starting in | f ) <8 |0) (in- 
stead of in \$i n ) = |J.)8> |0) as in Eq. {SJ|) and iterating 
the walk in Eq. (|14|) more we see that the walk starting 
in I f) ® |0) has a drift to the right side, exactly opposite 
to the walk in Eq. I|16(l . To obtain a symmetric distri- 
bution we can start the walk in a superposition of | 4) 
and I I) and make sure that these two drifts do not inter- 
fere with each other. This can be achieved starting in the 
symmetric state \$ S ym) = T)+i| |))<8|0). Since the 
Hadamard walk does not introduce any complex ampli- 
tudes the trajectories from | |) will stay real and the ones 
from I I) will be purely imaginary; they will hence not 
interfere with each other, making the total distribution 
symmetric. 

Another solution to eliminate asymmetry of the walk 
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FIG. 6: The probability distribution obtained from a com- 
puter simulation of the Hadamard walk with a symmetric 
initial condition |l5|. The number of steps in the walk was 
taken to be 100. Only the probability at the even points is 
plotted, since the odd points have probability zero. 



is to use a different (balanced) coin, namely 



Y 




(17) 



It is not hard to see that this coin treats | |) and | J.) in 
the same way and does not bias the walk, independently 
of its initial coin-state. 

Figure shows the probability distribution on the po- 
sitions of a symmetric quantum walk. 

The pattern of the probability distribution is very intri- 
cate - a signature of the quantum world. The multitude of 
oscillations makes it hard to analyze the moments of this 
walk precisely. Using two different approaches - combina- 
torial techniques (recursions, path counting - initiated by 
Meyer 0) - and a more physical path-integral approach 
(which expresses the walk in the Fourier domain in terms 
of integrals that are amenable to asymptotic analysis) 
Ambainis et al. |l2j| give an asymptotic analysis of the 
variance of the quantum random walk|56fl. 

As mentioned before the classical symmetric random 
walk on the line after T steps has a variance a 2 — T, so 
the expected distance from the origin is of order a = \/T . 
By contrast it can be shown that the quantum random 
walk has a variance that scales with a 2 ~ T 2 , which 
implies that the expected distance from the origin is of 
order a ~ T - the quantum walk propagates quadratically 
faster! 

Furthermore the walk spreads roughly uniformly over 
the positions in the interval [— ^= , - j=] as can be seen 

from Fig. and shown analytically 12] . This is again 
in stark contrast to the classical case in which the dis- 
tribution is peaked around the origin and drops off ex- 
ponentially several standard deviations a away from the 
origin. 




FIG. 7: The classical symmetric random walk on the line with 
one absorbing boundary placed in the origin. The starting 
point of the walk is in position 1. 



To uncover more striking differences of the quantum 
walk we can place it on a bounded line, either one-sided 
or two-sided |l2l |l6( . In other words we can insert one or 
two absorbing boundaries on the line. Formally an ab- 
sorbing boundary in position |6) corresponds to a partial 
measurement of the process at every time step. More pre- 
cisely the unitary step U of the walk will be followed by a 
measurement Mb on the position space Tip described by 
the two projections onto \b) and B± (the space orthogo- 
nal to \b)): 



MM) = 



\b) 

-jbW\b) 



Pb -- 

PB A 



\m\ 2 
= i-\m\ 2 



(18) 



where pb is the absorption probability and pb±_ its com- 
plement. For example if the current state of the system 
before measurement is given by 



I*) 



1 



(2|t>®|0)-|t)®|l)+3|i>®|l» (19) 



then after the measurement Mq the state will be | t) <8> |0) 
with probability 2/7 (in which case we say absorption 
occurred and stop the walk) or otherwise in the state 
^=(-| T) <8> |1> + 3| |) ® |1» with probability 5/7. 

Now define one step of the random walk with absorbing 
boundaries as U (Eq. H15l) ~) followed by Mb (or Mb and 
M^ in the case of two absorbing boundaries placed in |&) 
and | b')). If the measurement result of Mb (resp. My) 
gives \b) (resp. |&')) the walk is stopped, otherwise the 
next iteration is applied. 

Let us review what is known about a classical random 
walk with one absorbing boundary. For concreteness let 
us assume that the walk is started in position 1 and that 
the boundary is placed in position b = 0, as in Fig. \7\ 

It is well known that the probability p to ever get ab- 
sorbed by the wall in the origin is p = 1. This is very 
easy to see via a recursive reasoning on p: starting in 1 
the walk hits with probability 1/2, in which case it gets 
absorbed. Otherwise (with probability 1/2) its position 
is 2. The probability P20 to ever hit from 2 is the prob- 
ability P21 to ever hit 1 from 2 times the probability pio 
to ever hit from 1. Both P21 and pio are equal to p 
(the walk is homogeneous in space), which leads to the 
recursion 



11 2n 
P= 2 + 2 P2lPw = 2' P ' 



(20) 
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with the solution p — 1. In the classical world a ran- 
dom walk always gets absorbed by a boundary and never 
escapes to infinity 

Not so for the quantum walk. Careful analysis [l2l ITr| 
shows that the quantum Hadamard walk with absorbing 
boundary in |0), starting in the state | f) <8> |1) has a 
non-vanishing probability to escape. The probability to 
be absorbed at |0) turns out to be p qU antum — 2/7r 
0.63662, which gives an escape probability of w 0.36338. 
This behavior is general: for different initial states and 
distance to the boundary the numerical value for the ab- 
sorption probability p qU antum changes, but the escape 
probability is non-zero and finite in all cases. 

This is only the tip of the iceberg. After the semi- 
nal papers by Ambainis et al. 0] and Aharonov et al. 

a lot of research has been focused on quantum ran- 
dom walks uncovering more and more results of quantum 
'weirdness'. We will survey more of them in the following 
Section. 

In the spirit of the random walk on the line we can gen- 
eralize the definition to quantum random walks on gen- 
eral graphs. Various aspects of quantum random walks 
on graphs and in hig her dimensions have recently been 
studied in [H[ll El S El E3 To set the stage let us 
define the random walk for d-regular graphs first. These 
are graphs with vertex-degree d, i.e. each vertex has d 
outgoing edges. The coin Hilbert space He is of dimen- 
sion d. For every vertex let us internally label all the 
outgoing edges with j e 1 . . . d such that each edge has 
a distinct label. Let us call e° v an edge e = (v, w) which 
on u's end is labeled by j. The state associated with a 
vertex v pointing along an edge labeled j is | j) <g> \v) (cor- 
responding to | |) <£> |i) on the line or circle). Now we can 
define a conditional shift operation S 



S\j) ® \v) 



\j) <X> \w) ifej = (v,w) 
o/w 



(21) 



S moves the particle from v to w if the edge (v, w) is 
labeled by j on u's side. If the graph is not d-regular 
we can still define S in the same way if we let d be the 
maximum degree of any vertex in the graph. Figure |S] 
shows a graph and a valid labeling. 

The coin-flip C is now a d-dimensional unitary trans- 
formation, which gives a lot of freedom in the choice of 
the coin-flip operation. The choice of C will have to be 
guided by the specifics and symmetries of the random 
walk we wish to obtain. We have already seen that even 
in the simple case of a walk on the line the choice of the 
coin operator can have a serious impact on the shape of 
the resulting distribution, and this effect is only amplified 
in higher dimensions (see 0] for examples of the effects 
of various coins on a two- and three-dimensional grid) or 
on graphs of higher valence (degree) . 

If the graph is not regular, i.e. if the vertices have 
varying degrees, we can use one of several tricks to define 
a random walk with a coin. Either we add self-loops to 
each vertex of degree less than d (d is the maximal degree) 




FIG. 8: A graph with various degrees and a labeling of the 
edges for each vertex. Note that the label associated with and 
edge (v, w) can be different from d's end and from w's end. 
We have added self-loops to the graph to make it regular so 
that a single 4-dimensional coin can be used. 



and include them into the labeling (with the same label 
on both its ends), as in Fig. |SJ The shift applied to a 
self-loop will just keep the walk in place. Another option 
is to keep the irregularity of the graph and use different 
coin operators Cd' (of different dimensions) for vertices 
of degree d' < d. In that case the coin-operation has to 
be conditioned on the position of the particle and can no 
longer be written in the separable form of Eq. I|15|) 58] as 
(C<8)/). However whenever we can define the correspond- 
ing classical walk we are able to define the quantum walk 
and implement it efficiently on a quantum computer (see 
Sec. ITvTl . 

If we wish to retain the property in general graphs that 
the coin is balanced (i.e. every direction is obtained with 
equal probability if we measure the coin space) we can 
use the following coin which generalizes H 



/I 



DFT = 73 



1 



1 



1 \ 



\1 



ud-l w 2(d-l) 



,d-l 



(22) 

Here ui = exp(27ri/d) is a d-th root of unity. Clearly 
the unitary DF T-coin|59j transforms each direction into 
an equally weighted superposition of directions such that 
after measurement each of them is equally likely to be 
obtained (with probability 1/d). 

Let us give an example where a non-balanced coin is 
used. The graph is the hypercube of dimension d. Figure 
[5] shows the hypercube in 3 dimensions. The vertices of 
the d-dimensional hypercube can be enumerated by li- 
bit strings of 0's and l's. Two vertices given by two bit- 
strings arc connected whenever they differ by exactly one 
bit (e.g. 001101 and 011101). If we define the Hamming- 
distance du{x,y) between two bit-strings x and y to be 
the minimum number of bits that need to be flipped to 
obtain y from x; then the hypercube connects all vertices 
of Hamming distance dn = 1- The number of l's in a 
bit-string is called its Hamming weight. 

Note that the labeling of the edges of the hypercube 
can be chosen such that for each edge the labels from 
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FIG. 9: The hypercube in d — 3 dimensions. Vertices cor- 
respond to 3-bit strings. Edges are labeled by 1,2,3 (boxed) 
according to which bit needs to be nipped to get from one 
vertex of the edge to the other. 



both sides coincide 60] . 

The classical simple random walk on the hypercube 
(where in each step a neighboring vertex is chosen with 
probability 1/d) has a high symmetry. If we start the 
walk in the vertex 00 ... 0, for example, then all d vertices 
of Hamming weight 1 can be interchanged without mod- 
ifying the random walk. Similarly all vertices of equal 
distance to the starting vertex, i.e. of equal Hamming 
weight, can be interchanged without changing the walk. 
This implies that all vertices of same Hamming weight 
have the same weight in the probability distribution of 
the random walk. This allows us to 'cumulate' all the ver- 
tices of same Hamming weight into a single vertex and to 
reduce the symmetric random walk on the hypercube to 
a biased random walk on the line, as is shown in Figure 
1101 for the hypercube of dimension 3. Such a reduction 
to the line will substantially simplify the analysis of the 
random walk when we are interested in questions like the 
speed of propagation of the random walk to his opposite 
corner. 

If we wish to retain this symmetry property for the 
quantum random walk it can be shown |17| that the only 
coins that are unitary and permutation symmetric are of 
the following form 



G a ,b — 
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b . 


. b 
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b . 
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. b 
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b . 
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\b 
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b . 


. b 


a) 



(23) 



where a and b are real, 1 — § < \a\ < 1 and b = ±(1 — a). 
Among all these coins we will pick the coin G with 
a = 4 — 1 and b= \ and will sometimes call it 'Grover- 
coin'|61J . Among all coins G a ,b this is the coin which is 
the farthest away from the identity operator (Gi,o = I). 
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FIG. 10: The classical simple random walk on the 3- 
dimensional hypercube reduced to the walk on a line of 4 
points. In the original walk the transition probabilities are 
1/3. If we bundle nodes of same Hamming weight wh to a 
single point on the line we obtain a biased random walk on 
the line. The new transition probabilities for the walk on 
the line are indicated. For the walk on the d-dimensional hy- 
percube we can reduce the graph to a line of d + 1 vertices. 
The transition probability from vertex i to vertex i + 1 in the 
general case is given by = 2=i and from i + 1 to i by 



G is not a balanced coin because the probability to not 
change directions (p — (1 — |) 2 ) is different from the prob- 
ability to flip to one of the d— 1 other directions [p = -4). 
If we start the quantum walk in an equal superposition 
of all the directions and measure the coin space TLc after 
each iteration the resulting (classical) walk will have a 
higher propensity to go back and forth on the same edge 
than to switch directions. However the permutation in- 
variance will allow us to reduce the quantum walk on the 
hypercube to a quantum walk on the line. This would 
not be possible if we used the d-dimensional DFT-coin, 
for example. We will return to this graph in Sec. II VI 



B. Continuous time quantum random walk 

The discrete time model of quantum random walks is 
only one way to intersperse quantum effects with random 
walks. Another route to generalize the phenomenon de- 
scribed in Section [H] has been taken by Farhi et al. in 
1998 0jE3- At the outset this approach seems very far 
from the discrete time random walk we have described 
above. Later, similarities will appear, however. 

The continuous time random walk takes place entirely 
in the position space Tip. No coin space is needed and no 
coin is flipped. The intuition behind this model comes 
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from continuous time classical Markov chains. Let us 
illustrate this with the simple classical random walk on a 
graph with vertex set V. A step in the classical random 
walk can be described by a matrix M which transforms 
the probability distribution over V. The entries M^j give 
the probability to go from i to j in one step of the walk. 
Let p 1 = (p\,p t 2 , ■ . ■ be the probability distribution 

over the vertices of V at time T. Then 



E 



which means 



Pt+i 



(24) 



(25) 



The entries of are non-zero only if there is a non- 
zero probability to go from i to j, i.e. when i and j are 
connected. The entry Mij is the probability to go from i 
to j and is equal to -j- (in the so called 'simple' random 
walk) where di is the degree of i. For instance the matrix 
M for a simple random walk on a circle of N vertices is 
given by 



M = 



(o 1 
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(26) 



...,0), cor- 
responding to a starting position of 0, to p 1 = 
(0,|,...,0,|), p 2 = (i,0,|,0,...,0,i,0), f = 
(0, |, 0, |, 0, . . . , 0, g, 0, |) etc. corresponding to the dis- 
tributions as listed in Table |3J 

The process given by the iterations of M transforms 
the state at integer times only. To make the process con- 
tinuous in time we can assume that transitions can occur 
at all times and the jumping rate from a vertex to its 
neighbor is given by 7, a fixed, time-independent con- 
stant. This means that transitions between neighboring 
nodes occur with probability 7 per unit time. The in- 
finitesimal generator matrix H of such a process is given 
by 

-7 i j and i and j connected 
= { i 7^ j and i and j not connected (27) 
da i = j 

If Pi (t) denotes the probability to be at vertex i at time 
t then the transitions in analogy to Eq. (|24|) can be 
described by a differential equation 



dpijt) 
dt 



(28) 




FIG. 11: The graph G4. Two binary trees of n = 4 levels are 
glued together at their leafs in a symmetric fashion |2Rl |. 



Solving the equation we obtain p(t) = cxp(— Ht)p(0) in 
analogy to Eq. (|25|l . Note the similar structure of H and 
M . In the theory of classical Markov chains many con- 
nections between the discrete and continuous time model 
can be made and many of the quantities of interest, like 
mixing times (see Sec. lIVBf) and ab sorp tion probability, 
have similar behavior in both cases l2fl . 

The idea of Farhi and Gutmann |24j was to carry this 
construction over to the quantum case. Their key idea is 
that the generator matrix will become the Hamiltonian 
of the process generating an evolution U(t) as 



U(t) = exp(-iHt). 



(29) 



If we start in some initial state l^m), evolve it under 
U for a time T and measure the positions of the result- 
ing state we obtain a probability distribution over the 
vertices of the graph as before. 

With this definitions in place Farhi and Gutmann |24| 
study the penetration of graphs by the quantum random 
walk. In |25j a finite graph is given where classical and 
quantum walks give an exponential separation in expected 
hitting time 62] . We will give this example here to illus- 
trate yet another aspect of the difference of classical and 
quantum random walks, this time regarding the expected 
hitting time of certain nodes in a graph. The expected 
hitting time between two points S and T in a graph is 
the time it takes the random walk on average to reach T 
starting at S. The graph G n consists of two n-level bi- 
nary trees glued together at their leafs as seen in Figure 

EH 

We start the random walk at the root of the left tree. 
Rather than the hitting time we study a slightly modified 
related quantity: the probability of hitting the root of the 
other (right) tree as a function of time. In other words, 
we are interested in how long it takes to propagate from 
the leftmost vertex to the rightmost vertex as a function 
of n. 
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Consider the classical random walk first. The ver- 
tices of G n can be grouped in columns indexed by j £ 
{0, 1, . . . , 2n}. Column contains the root of the left 
tree, column 1 contains the two vertices connected to 
that root, etc. Note that column n contains the 2™ ver- 
tices in the middle of the graph and column 2n is the 
root at the right. 

As before in the example of the hypercube (Fig. I10|) 
the symmetry of the walk allows us to reduce it to a 
random walk on a line. We need only keep track of the 
probabilities of being in the columns. In the left tree, 
for < j < n, the probability of stepping from column 
j to column j + 1 is twice as great as the probability of 
stepping from column j to column j — I. However, in the 
right tree, for n < j < 2n, the probability of stepping 
from column j to column j + 1 is half as great as the 
probability of stepping from column j to column j — 1. 
This means that if you start at the left root, you quickly 
move to the middle of the graph, but then it takes a 
time exponential in n to reach your destination. More 
precisely, starting in column 0, the probability of being 
in column 2n after any number of steps is less than 2~ n . 
This implies that the probability of reaching column 2n 
in a time that is polynomial in n must be exponentially 
small as a function of n. 

We now analyze the quantum walk on G n p5l | starting 
in the state corresponding to the left root and evolv- 
ing with the Hamiltonian given by (PZJ. With this ini- 
tial state, the symmetries of H keep the evolution in a 
(2n + l)-dimensional subspace of the (2™ +1 + 2" - 2)- 
dimensional Hilbert space. This subspace is spanned by 
states \j) (where < j < 2n), the uniform superposition 
over all vertices in column j, that is, 



1 



3) = 



E 



a e column j 



where 



J 1 2 2n -i n<j<2n. 



(30) 



(31) 



In this basis, the non-zero matrix elements of H are 

(j\H\j±l) = -V2~f (32) 
/ 2 7 j = 0, n, 2n 



(3\H\j) 



37 otherwise, 



(33) 



which is depicted in Figure 1121 (for n — 4) as a quantum 
random walk on a line with 2n + 1 vertices. 

As a first attempt to solve this reduced problem we can 
approximate the walk on this finite line by a walk on an 
infinite and furthermore homogeneous line by extending 
the Hamiltonian and replacing all diagonal entries of 2 
by 3 as seen in Figure ED 

Solving this modified problem is the lore of a quantum 
physicist and the solution can readily be given in terms 
of Bessel functions (see [2tJ for a much more complete 
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FIG. 12: The elements of the Hamiltonian of the quantum 
random walk on G n when reduced to the line for n = 4 |25|. 
For convenience 7 is set to 1. 



FIG. 13: Approximation of the finite line by an infinite ho- 
mogeneous line 12511 . The diagonal elements of H are all set 
equal to 3. The off diagonal elements connecting the vertices 
arc still — v/2. 



treatment). It can be seen that the speed of propagation 
of the random walk on the infinite homogeneous line is 
linear in the time T . In other words there are constants 
a < b such that when T € [an, bn] then the probability 
to measure the walk a distance 2n + 1 from the staring 
point (i.e. at the root of the right tree starting from the 
root of the left tree) is of order 1/y/n. 

If we account for the finiteness of the line we are facing 
the standard quantum mechanical problem of a particle 
in a square-well potential, which is still not hard to solve. 
It is not entirely obvious that the finiteness of the line 
and the 'impurity' at level n do not change this behavior 
significantly Numerical studies (2^ confirm that the 
linear propagation behavior is not significantly changed. 
To get a hand on the exact hitting time of this walk very 
recently |27| an upper bound for the propagation time is 

found. The precise statement is that if a time T is chosen 

4 

randomly T S [0, ^-j and the walk is run for a time T, 
then the probability to measure the state in the root of 
the right tree is phu > r(l — e). For algorithmic purposes 
we can repeat the random walk a polynomial number of 
times to almost surely find the particle at the right root 
at least once. 

We can also implement a discrete time random walk on 
the graph G n using the 3-dimensional Grover coin G (Eq. 
1)23(1) at the vertices of degree 3 and the 2-dimensional 
one at vertices of degree 2 (wh ich just flips the two direc- 
tions). Numerical studies [23 indicate that this discrete 
time random walk shows an extremely similar behavior 
to its continuous counterpart. In particular it is possible 
to reduce the walk to a discrete time walk on the line on 
which the walk propagates in linear time. This is just 
one example of the similarity of both models, of which 
there are many others. 

Before we give more results on the difference in behav- 
ior of classical and quantum random walks let us pause 
for a moment to investigate the importance of random 
walks for algorithms in computer science. 
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IV. RANDOM WALKS IN COMPUTER 
SCIENCE 

Classical random walks form one of the cornerstones 
of theoretical computer science 1^3, EI ■ As algorithmic 
tools, they have been applied to a variety of central prob- 
lems I32L such as estimation of the volume of a convex 
body |33], approximation of the permanent of a matrix 
[34|. and discovery of satisfying assignments for Boolean 
formulae |35l | . They provide a general paradigm for sam- 
pling and exploring an exponentially large set of com- 
binatorial structures (such as matchings in a graph), by 
using a sequence of simple, local transitions |32j ] . 

The behavior of many an algorithm that use random 
walks depends on quantities like its mixing time (the time 
it takes the walk to approach its stationary distribution; 
this is important if we wish to sample from a set of com- 
binatorial structures according to some distribution cf. 
the Markov Chain Monte Carlo Method 64] ) or its ex- 
pected hitting times between two vertices of the under- 
lying graph. We will illustrate the latter point with two 
simple examples before we give the more formal defini- 
tions. 



A. Two examples 

1. S-T Connectivity 

Let us illustrate how random walks are harnessed in 
classical computer science by describing how they help to 
solve the s — t connectivity problem. Given an undirected 
graph G with vertex set V of size |V|, the problem is 
to decide whether two vertices s and t are connected or 
not. This problem is a natural abstraction of a number 
of graph search problems and is important in the study 
of space bounded complexity classes in computer science 

m 

It is easy to see that a standard graph search algorithm 
that keeps track of all the vertices the search has visited 
can do the job in a time linear in the number of edges 
of the graph. Let us assume now that our computer 
is limited in memory and can keep track only of a few 
vertices at a time|65|. It can perform a random walk on 
the vertices of the graph starting in s. At each step of 
the walk the computer choses one of the adjacent edges 
with equal probability and moves to the vertex that is 
connected to the current position by that edge. If the 
computer reaches t in the course of the walk it outputs 
YES. If it never sees t after T = \V\ 3 steps it outputs 
NO. 

Clearly the machine will never output YES when s and 
t are not in the same connected component. It turns out 
that the simple random walk on a graph with \V\ ver- 
tices has at least probability 1/2 to hit any given vertex 
starting from any other in \V\ 3 steps. This means that 
our machine will err with probability 1/2 when it out- 
puts NO. We can make this probability arbitrarily small 



by repeating the algorithm several times. Assume that 
s and t are connected. Then the error probability of an 
algorithm that iterates the above procedure k times is 
simply the probability to output NO k times. Since we 
assume that each iteration is independent of the previ- 
ous one, this probability is given by (l/2) fe . If we want 
to make this probability very small (e) we need to repeat 
the algorithm roughly k = log(i) times. 

2. 2-SAT 

Let us give another computational problem that can 
be solved by random walks on graphs, 2-SAT. Let us 
define more generally the class of decision problems called 
SAT. In an instance of SAT we are given a set of logical 
clauses C±, . . . C m . The Boolean inputs are the variables 
X1.X2, ■ ■ ■ X n which can take the values or 1 . They can 
appear in either uncomplemented (Xi) or complemented 
(-iAj) form in the clauses. A clause is said to be satisfied 
if at least one of the variables in it is true (an unnegated 
variable Xi is true if Xi = 0, a negated variable ~^Xi is 
true if Xi — 1). For example the clause C = X\ V ~^X 2 V 
A 3 is not satisfied only if XiX 2 X 3 = 101. We want to 
know if all the clauses can be satisfied simultaneously by 
some assignment of and 1 for X\ , . . . X n (then we say 
that the formula C\ A C2 A . . . A C m is satisfiable and 
output YES). The class 2-SAT is the class of all SAT 
decision problems where each clause contains only two 
variables. For example take the following instance of 2- 
SAT on X 1 ,X 2 ,X 3 : 

$(X 1 ,X 2 ,X 3 ) = (Xi V -1X2) A (->Xi V A3) 

A(X 2 V X 3 ) A (-.JTi V ^X 3 ) (34) 

The underlying decision problem has the answer YES 
and a satisfying assignment is X1X2X3 = 110. 

The class of satisfiability problems is very important 
in theoretical computer science [30| . For instance in a fa- 
mous piece of work (the Cook-Levin theorem) it has been 
shown that 3-SAT is as hard as any problem in the class 
NP. This class contains all problems that can be solved in 
so called non-deterministic polynomial time. That means 
that if we guess a solution then we can verify that it is 
indeed a solution in polynomial time. For instance if we 
were presented with a guess for the satisfying assignment 
of a formula in 3-SAT we are able to verify this with a ma- 
chine running in time polynomial in the number of vari- 
ables n by checking all the clauses Ci, . . . C m . A plethora 
of decision and optimization problems are in the class NP. 
The Cook-Levin theorem states that if we had an efficient 
algorithm to solve 3-SAT (i.e. to find an answer YES/NO 
to the decision problem) we can change this into an algo- 
rithm to solve any other problem in NP with only poly- 
nomial overhead in time. Among the variety of problems 
known to be in NP we can find Shortest Common Su- 
perstring (relevant for gene-reconstruction from sequenc- 
ing in biology) , the Traveling Salesman Problem (impor- 
tant for tour-optimization), Bin Packing problems, the 
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FIG. 14: A random walk on assignments to the 2-SAT formula 
of Eq. (I34t . Each step increases or decreases the Hamming 
distance to the true assignment 110 by 1. This gives rise to a 
random walk on the line. 

Minesweeper game and Protein Folding, to name just a 
few. It is not known if a polynomial time (in n) algorithm 

to solve 3-SAT exists (the famous P = NP problem). 

Let us now give a random walk algorithm that solves 
2-SAT in time proportional to n 2 . To simplify things 
imagine that our formula is satisfiable and that it has 
a single truth assignment (in the example of Eq. l|3"4"|) 
this is 110). The algorithm picks a random assignment 
of 0, 1 for the variables X\, ■ ■ ■ , X n . It checks the clauses 
one by one if they are satisfied until it either runs out 
of clauses (this means the assignment was satisfying, the 
algorithm stops and outputs YES) or it finds a clause C 
that is not satisfied. Say the two variables in the clause 
are Xi and Xj . In the true assignment one of these two 
variables must be flipped (0 <-> 1). Now the algorithm 
choses one of the two variables Xi or Xj randomly and 
flips it. This is the new assignment and the algorithm 
restarts anew checking the clauses one by one etc. Let 
us look what happens to the Hamming distance between 
the true assignment (110 in the example of Eq. (JH2J) an d 
the current assignment of the algorithm when we come to 
flipping Xi or Xj . In the true assignment one of the vari- 
ables Xi, Xj must be flipped (since it satisfies the clause 
C). So with probability 1/2 the algorithm flips the right 
variable. In this case the Hamming distance between the 
new assignment and the true assignment has decreased 
by 1. Or the algorithm flips the wrong variable, in which 
case the Hamming distance increases by l|6g. So we 
are faced with a random walk on the Hamming distances 
with transition probabilities 1/2. It is well known that 
such a walk on the line of length n will hit the corner with 
high probability after T ~ n 2 steps. Figure ITU illustrates 
the algorithm for the example of Eq. i|34[l • 



Although this random walk algorithm is not the op- 
timal one to solve 2-SAT it is very instructive. Many 
randomized algorithms can be viewed as random walks. 
The efficiency of these random algorithms is then closely 
related to the expected hitting time from a vertex (the 
initial state) to another (the solution). In fact the best 
known algorithm for 3-SAT [35[ relies on a sophisticated 
version of this random walk idea. 

If a quantum random walk could improve the hitting 
time on certain graphs this could give rise to efficient 
quantum algorithms for problems that cannot be solved 
efficiently on a classical computer. It is yet another in- 
centive to study the behavior of quantum random walks. 

B. Classical Random Walks 

Let us state without proofs some known facts about 
random walks on graphs |36j . 

A simple random walk on an undirected graph G(V, E), 
is described by repeated applications of a stochastic ma- 
trix M, where Afjj = -r if (i,j) is an edge in G and di 
the degree of i (we use the notation introduced in Sec. 
IIII Bp . If G is connected and non-bipartite, then the dis- 
tribution of the random walk, p 5 " = M T fP converges to 
a stationary distribution ir which is independent of the 
initial distribution jr. For a d— regular G (all nodes have 
the same degree d), the limiting probability distribution 
is uniform over the nodes of the graph. 

In many computational problems the solution can be 
found if we are able to sample from a set of objects ac- 
cording to some distribution 67] . This problem is often 
approached by setting up a random walk on a graph 
whose nodes are the objects we want to sample from. 
The graph and the walk are set up in such a way that the 
limiting distribution is exactly the distribution we wish 
to sample from. To sample we now start the walk in 
some random initial point and let it evolve. This type of 
algorithm is only efficient if the random walk approaches 
the limiting distribution fast. 

There are many definitions which capture the rate of 
convergence to the limiting distribution. A survey can be 
found in ^5% . A frequently used quantity is the mixing 
time given by 

M e = mm{T\Vt>T,f : - vr 1 1 < e}, (35) 

where we use the total variation distance to measure the 
distance between two distributions p, q: 

Il#-<z1l=£^-e (36) 

i 

It turns out that the mixing time is related to the gap 
between the largest eigenvalue Ai = 1 of the stochastic 
matrix M, and the second largest eigenvalue A2 in the 
following way |37l ] 



(l-A 2 )log2e " 
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< 71 TT ( max lo § ^ 1 + lo S e 1 537) 

(1 — M) 1 

This powerful theorem provides a very useful connection 
between mixing times and the second largest eigenvalue 
A2 of the transition matrix M. In particular it shows 
that A2 is the only eigenvalue of M that really matters 
to determine the mixing behavior of the walk. This will 
be very different in the quantum case. 

Circle: It is well known that for the simple random 
walk on an N— circle, the mixing time is quadratic, M e ~ 
N 2 ■ log(l/e). Similarly the expected time T we need to 
run the walk such that starting at a vertex i we hit a 
certain vertex j with probability close to 1 is T ~ N 2 . 

Hypercube: For the d-dimensional hypercube defined 
in Sec. IIII Al the mixing time scales with d as M e ~ 
dlogd log(l/e). The expected hitting time T from one 
corner of the hypercube, say 00 ... 0, to its opposite, 
11 ... 1, scales exponentially with the dimension d as 
T ~ 2 d . We will see a sharp quantum improvement here. 

C. Quantum Computers and Circuits 

We have developed an intuition for why quantum ran- 
dom walks could be useful for algorithms. Before we de- 
scribe some more results on the mixing and hitting times 
of quantum random walks let us first expand a little on 
the machine that would implement the quantum random 
walk - the quantum computer. We will only give a rushed 
collection of facts to acquaint the reader with the basics, 
for a very detailed introduction see 0. 

A quantum computer functions according to the laws 
of quantum mechanics. In particular its 'computations' 
have to be unitary transformations on the state space. 

Qubits: A classical computer operates on strings of 
bits. Similarly a quantum computer acts on qubits. Each 
qubit is a two level system spanned by the two states 
{|0), |1)} (or {| |), I J,)} to make a closer connection to 
spins). We can think of the state Hilbert space of n 
qubits to be the tensor product of n 2-dimensional qubit 
spaces Hf n . Another way to think of the states is to say 
that they are spanned by n-bit strings \x1X2 ■ ■ ■ x n ) with 
Xi G {0, 1}. A state of the quantum computer can be in 
any superposition of these basis states. For a physicist 
it might be easier to think about a collection of n spin-i 
particles. 

Gates: The computation of a classical computer can 
be described by a circuit consisting of gates (like AND or 
NOT) that act on one or two bits. The computer is ini- 
tialized in some input state and computes some Boolean 
function of it. Similarly a quantum computer can be 
modeled by a circuit consisting of unitary gates that act 
on one or two qubits only. The goal is to compute some 
unitary transformation on the initial state. 

It is well known that every Boolean function can be de- 
composed into a succession of gates from a universal set 
(like AND and NOT). In other words any computation 
on a classical computer can be written as a succession 



of AND and NOT. One of the achievements of quantum 
computing theory has been the insight that any unitary 
computation on qubits can be decomposed into a succes- 
sion of quantum gates that act on one or two qubits only. 
One such universal quantum gate set is given by single 
qubit unitary matrices U S i ng i e (these are 2 by 2 matri- 
ces giving all the generalized rotations of a single spin-i 
system) together with the so called CNOT 
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To build a quantum computer then means to build a de- 
vice that can implement single qubit unitaries together 
with a two-body transformation that implements the 

cnotIH. 

The question is whether such a quantum computer 
with its unitaries can even perform the computations of 
a classical computer. After all, unitary transformations 
are reversible, so all the quantum computations are re- 
versible and the class of these reversible computations 
might not contain some classical circuits. It has been 
shown by Bennett [38j that this is not a serious problem: 
Any classical circuit can be made reversible by replacing 
its gates by a reversible three bit gate and adding some 
more scratch space. Furthermore these reversible gates 
can be implemented unitarily by a quantum computer. 
In particular this implies that a quantum computer can 
emulate any classical computation efficiently: quantum 
computers are at least as strong as their classical broth- 
ers. 

The added bonus of using the laws of quantum me- 
chanics is that the unitary power might be strictly larger 
than the classical one. The quest for quantum algorithms 
is exactly to explore this added power. 

Returning to random walks we can see that a quan- 
tum computer can efficiently implement any discrete time 
quantum random walk whenever a classical computer can 
implement the underlying classical random walk. This 
is because we can take the classical circuit that imple- 
ments the classical random walk, turn it into a quan- 
tum circuit (after making it reversible) and replacing the 
classical coin-flip by a quantum coin-flip, which can be 
done efficiently on the small space of coin-qubits. That 
means that if we take a classical random walk at the basis 
of some classical algorithm and show that the quantum 
random walk improves its behavior, we will immediately 
have a quantum algorithm with some speed-up over the 
classical. 

The situation is not quite as easy for the model of con- 
tinuous time quantum random walks. It is not a priori 
clear how to adapt its continuous evolution to the discrete 
quantum circuit model. However it has been shown re- 
cently mm that simulation of continuous time random 
walks can be done efficiently on a quantum computer. 

With this in place we are now set to explore some more 
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differences between classical and quantum random walks. hence 
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D. Results from Quantum Random Walks 

1. The circle 

The first and maybe easiest finite graph to explore with 
quantum random walks is the circle of N vertices |llj . 
This example has most of the features of walks on general 
graphs. 

Recall that any classical random walk approaches a 
stationary distribution independent of its initial state 
the classical random walk looses its memory 

All 



(Sec. HVBt 

This cannot be true for the quantum random walk, 
transformations are unitary and hence reversible, there- 
fore the walk can never lose its recollection of the initial 
state and therefore it does not converge to a stationary 
distribution. To be able to speak of mixing toward a sta- 
tionary distribution we have to introduce some 'forget- 
ting' into the definitions. There is a standard way to do 
this, the so called Cesaro limit. Instead of studying the 
probability distribution jf induced by the random walk 
after a measurement at time t and its limit as t grows we 
pick a time s, < s < t, uniformly at random, i.e. our 
probability distribution c* is an average distribution over 
all measurement results between 1 and i: 



1 * 
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With these definitions in place it is not very hard to see 
that c* converges to a stationary distribution. We will 
show this below for the Hadamard walk on the circle 
[Tl| - also to give an idea about the kind of calculations 
involved in these types of problems - but the result is 
true for all finite graphs. 

In order to analyze the behavior of the quantum walk, 
we follow its wave-like patterns using the eigenvectors 
and eigenvalues of the unitary evolution Ut- To ob- 
serve its classical behavior, we collapse the wave vec- 
tor at time t, |* t ) — J7*|*q), into a probability vector 
fr . The probability to measure the particle in position i 
at time t, $ = |(t,i|* t >| 2 + |(|,i|* t )| 2 , where we write 
| short for | |) ® \i) and similarly for | i,i). Let 

{(Afe, \vk)) : k = 1 . . . 2N} be the eigenvalues and eigen- 
vectors of U. Note that \* k = A^T 1 because U is unitary. 
Expanding the initial state |^o) — Y^t=i a k\ v k) we ob- 
tain |* t ) = 17*1*0) = J2l=i a kK\v k )- Putting this all 
together we get for the ith component 
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In the limit t — > oo we have 
t 



Afe — A; 
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t(l-(A h A*)) _ Afc ^ A; 



(41) 



-^oo ^ ^ a k ai(a,i\v k )(vi\a,i) =: n l . (42) 



a=U k,l = l 

Afc = A; 

In the case of non-degenerate eigenvalues Afc, Eq. H42|) 
simplifies further to Wi = Yla=1,l Sfc=i \ a k\ 2 \(a,i\v k )\ 2 ■ 

Eqs. H4U|) and (|42[l already reveal some crucial differ- 
ences between classical and quantum walks. Whereas in 
the classical case the stationary distribution jfi is inde- 
pendent of the initial state of the walk, in the quantum 
walk this is not the case in general, as is manifested by the 
presence of the expansion coefficients a k in the expres- 
sion for TTi. However in the case of the circle (and many 
other graphs from a family called Cayley graphs [63) the 
stationary distribution is uniform and independent from 
the staring state as long as this state is localized in one 
position |11| . 

Another remarkable difference is the dependence of the 
mixing time - the rate of convergence to 7? - on all eigen- 
values of U. Remember that in the classical case the 
mixing time was governed only by the second largest 
eigenvalue of the transition matrix M (see Sec. II V Bl . 
As can be seen from Eq. H4U|) the rate of convergence to 
7? in the quantum case is governed by the terms 
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All of these sums over pairs of eigenvalues enter Eq. I|40|l , 
and so the mixing time is determined by all | Afc — A; | . 

It is possible to bound the mixing time M e of the quan- 
tum random walk on the circle as M e < — l °f — with 
some more lines of algebra This gives a nearly 

quadratic speed-up over the classical walk on the circle, 
which mixes in time proportional to N 2 (see Sec. IIV B(l ! 
We have found yet another 'quantum advantage': any al- 
gorithm that uses a random walk on a circle can be made 
quadratically faster on a quantum computer! 

We might ask if there are graphs for which we can 
achieve even more drastic speed-ups. In computer sci- 
ence the transition between an easy and a hard problem 
is made between polynomial time and exponential time. 
Problems that take exponential time to solve are hard 
(the resources to solve them do not scale well with their 
size), problems that take only polynomial amount of time 
are 'easy' [23 ■ The interesting question is now: Can we 
find graphs such that the quantum walk mixes exponen- 
tially faster than its classical brother? Unfortunately the 
answer is 'NO' in all but possibly very contrived cases. 
As shown by Aharonov et al. Ill] the quantum walk 
mixes at most quadratically faster |7lj. 

Yet this is not the end of all hopes for exponential 
speed-up, as we will see in the next section. 



2. The hypercube 



It is possible to calculate the eigenspectrum of sev- 
eral matrices U corresponding to a random walk on a 
graph with some coin C. Moore and Russell [l7] anal- 
ysed the mixing time of the walk on the hypercube (see 
Sec. Illlf) along the lines of the previous paragraph, both 
in the discrete time case as well as in the continuous one. 
Their results for the mixing time are not very encour- 
aging, though. The hypercube is one of the graphs that 
mixes fast to the uniform distribution in the classical case 
already (it belongs to a class of so called expander graphs 
which all share this property), namely M e ~ dlogd (Sec. 

MM- 

The discrete time quantum walk turns out to do worse: 
Its mixing time M e is at least d 3 / 2 /e. Even more dramat- 
ically the continuous time quantum walk on the hyper- 
cube does not converge to the uniform distribution at 
all. 

We have mentioned another quantity of crucial interest 
in the study of random walks - the expected hitting time. 
We have shown the importance of fast hitting to solve 
some algorithmic problems in Sec. IIVI Furthermore it is 
here that we may hope to achieve great quantum advan- 
tage: the classical random walk on the hypercube takes 
exponential time to hit its opposite corner. 

In the quantum world the notion of hitting can be 
made precise in several ways (see for a systematic ap- 
proach). In order to know whether a walk hits a certain 
point measurements need to be made, yet we do not wish 
to measure the position of the walk completely to not kill 
its quantum coherences and make it classical. In [l^ two 
ways to define hitting time are presented: either we can 
wait for a certain time T (which we somehow determine 
in advance) and then measure the walk. If the probabil- 
ity to be in the node of interest is close to 1 we can call 
T a hitting time. To circumvent the problem of knowing 
T in advance we can alternatively use a definition similar 
to the absorbing wall described in Sec. IIHI At each step 
a measurement is performed which only determines if the 
walk is in the node of interest or not (Mb in Eq. (|18|l of 
Sec. IIII|I . If it is found there the walk is stopped. We will 
call T a hitting time if the walk has stopped before time 
T with high probability. In the classical case all these no- 
tions of hitting time from one corner of the hypercube to 
its opposite are exponential in d. Kempe jl8j shows that 
the quantum discrete time walk indeed behaves crucially 
different: The expected hitting time (both notions) from 
a corner to its opposite turns out to be polynomial in the 
dimension d\ This is an exponential separation of the 
classical and quantum behavior, similar to what we have 
seen on the graph G n for the continuous time quantum 
walk (Sec. IIII Bl) 72]. It rests to see how this speed-up 
can be harnessed algorithmically. 
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ENTRANCE 




FIG. 15: The modified graph of Fig. 1111 A big random cycle 
has been inserted between the two binary trees. This makes 
the classical task of penetrating the graph from root to root 
exponentially hard |27|l . 



3. An oracle separation 

The first successful attempt to make an algorithmic 
statement using exponentially faster quantum hitting 
times is given by Childs et al. ^fr. They provide a so 
called oracle separation using the example described in 
Sec. IIII Bl We have seen that the continuous time quan- 
tum walk penetrates the glued trees from root to root 
exponentially faster than the classical walk. This can 
be made into an exponential separation of classical and 
quantum query complexity in an oracle problem. This 
means the following: We are given an oracle (in form of 
a black box that we cannot open), which we can query 
with some input to get a specific answer. For instance 
in database search set-ups the oracle encodes a database 
with one marked item (like a telephone book where we 
happen to know one specific phone-number and wish to 
know who it belongs to). Given the input string x (a 
possible name) as a query, the oracle will reply if a; is 
not marked (doesn't belong to the phone number) or 1 
if it is the item we are looking for. Imagine that each 
such query costs a fixed amount; it is our task then to 
minimize the number of queries. It has been shown that 
a quantum computer (who can ask a query 'in superpo- 
sition') can find the marked item with only \f~N queries, 
where N is the number of items, whereas any classical 
algorithm will require of order N queries. 

Childs et al. construct an oracle which encodes a graph 
similar to the one in Fig. ^] It turns out that there 
are classical algorithms (random walks that keep track 
of the valence of vertices they visit) that can penetrate 
the graph of Fig. 1111 in polynomial time. To make the 
task classically hard the graph has to be modified by 
inserting a huge cycle between the trees, as in Fig. El 

The oracle now encodes the structure of the graph: Ev- 
ery node of the graph has a name unknown to us. We 
are given the name of the right root and wish to find the 



16 



name of the left root. Queried with the name of some 
vertex the oracle outputs the names of the neighbors of 
this vertex. Classically this oracle has to be queried ex- 
ponentially many times in n (number of levels of the tree) 
to obtain the name of the root of the left tree starting 
with the name of the root of the right tree. Intuitively, 
the best way to proceed classically is to ask the oracle 
for the names of the neighbors of the root. With the 
answer we ask for the names of the neighbors of one of 
these neighbors and so on making our way through the 
graph in the form of a random walk. We have seen that 
it takes exponential in n expected time for the classical 
random walk to hit the opposite root. Furthermore it 
can be shown |27| that this naive way of making our way 
through the graph is the best a classical algorithm can do. 
Using the continuous quantum walk|73fl to find the suc- 
cession of (quantum)-queries to the oracle the problem 
can be solved with only polyno mially many queries us- 
ing a quantum computer |27| |74| . For this (unfortunately 
not very practical) problem we have found an exponential 
separation between classical and quantum power based 
on quantum random walks 75j! 



4. What next? 

We hope that these first algorithms are only the begin- 
ning and that more algorithms can be found that employ 
the power of quantum random walks. For a recent search 
algorithm based on the discrete time walk on the hyper- 
cube see [22^. 

Many questions are still left open. We will mention 
some: For instance the quantum random walk on a cir- 
cle of even degree (and on many other graphs) does not 
converge to the uniform distribution (its eigenvalues are 
degenerate). Rather, its stationary distribution depends 
on the starting state. It would be interesting to know if 
one can harness this 'quantum memory' algorithmically. 

Another observation is that the exponential hitting 
time on both hypercube and glued trees seems to de- 
pend crucially on the symmetry of these graphs. Rapid 
hitting seems to indicate high symmetry and could be 
used to perform some sort of graph-'tomography', which 
might have algorithmic use. 

There are a lot of graphs that have not yet been stud- 
ied in the context of quantum walks (like the walk on 
the symmetric group, for instance, which is paramount 
in some interesting algorithmic questions like Graph- 
Isomorphism) . 

Furthermore the connection between the two quantum 
walk models is still not clear. In many cases both walks 
behave very similarly, in other cases there are differences. 
It would be useful to make the connection between these 
two models more precise. 



V. PHYSICAL IMPLEMENTATION 

We have attempted to convince the reader of the utility 
to study random walks - both for the interesting physical 
effects they show as well as for their potential algorith- 
mic use. In both cases it would be interesting to find a 
physical set-up to implement the quantum random walk. 
Of course a universal quantum computer can do that - 
as we have seen in Sec. IIVI However, this machine has 
yet to be built on a large scale. Random walks might 
exhibit specifics that could ease their implementation in 
certain physical systems, making use of the structure of 
the graph underlying the walk or the character of the 
coin, for example. Such an implementation does not nec- 
essarily constitute a full fledged quantum computer but 
could still be well suited to study random walks or to 
solve algorithmic problems based on them |22], |27j. In 
this spirit several proposals for various physical systems 
have been made E2, EH, • 

As of today it is not clear which physical architecture 
is the optimal one to build a quantum computer. As we 
have seen we need qubits and gates which can be imple- 
mented in a multitude of ways. This fungibility of quan- 
tum information gives rise to proposals ranging from nu- 
clear magnetic resonance systems, optical cavities, solid 
state architectures, optical lattices to ion traps, to name 
just a few. Similarly the options to implement a quan- 
tum random walk are vast. We will just touch upon a 
few, without giving a detailed description of the physical 
system; the reader is referred to the respective references 
for the details. 

Of course the positions of the walk \i) need not be real 
positions in the physical implementation: we can encode 
them into any discrete degree of freedom. Similarly the 
coin space need not correspond to the spin of a spin-i 
particle. 

For instance Travaglione and Milburn 42] suggest an 
ion-trap based implementation of the random walk on 
a line. In this set-up an ion is confined in a Radio- 
Frequency ion trap (see |45ll46| for an overview on the ion 
trap quantum computer). The positions \i) are encoded 
into the (discrete) motional states of the ion in the trap. 
The ion is laser-cooled to its motional ground state |0). 
The electronic internal states of the ion will encode the 
coin states | f) and | J,). Raman beam pulses are applied 
to implement a coin-flip. A displacement beam excites 
the motion correlated to the | f) internal state of the 
ion. Flipping the internal states with another Raman 
pulse and applying the displacement beam again effec- 
tively implements the conditional shift. This procedure 
can be iterated T times. To measure the positions of the 
walk after T steps it is possible to couple the motional 
state to the internal state of the ion and measure the lat- 
ter (see E2 f° r details). The main source of noise in this 
scheme is the degradation of the internal state of the ion 
(the coin state). 

With present day technology we might be able to im- 
plement a few steps of the random walk 76j . If nothing 
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else the variance of the quantum random walk so ob- 
tained can serve as a benchmarking protocol for ion trap 
quantum computers. 

Another recent suggestion by Sanders et al. 43] puts 
the random walk into a quantum electrodynamics device. 
The physical system is an optical cavity traversed by an 
atom. During the traversal the internal electronic levels 
of the atom interact with the cavity mode. Furthermore 
the internal atomic levels can be addressed and manipu- 
lated with a periodic sequence of microwave pulses. The 
focus of this proposal is to implement the random walk 
with existing present day technology in a microwave cav- 
ity. In this set-up the position of the random walk is en- 
coded into a single cavity mode, whereas the coin states 
are the internal atomic states [77j. Here a random walk 
on the circle is achieved. Let us sketch the scheme of 01 
in more detail (for the optical background see |47j'): 

The field mode of the cavity is described by the (in- 
finite) Hilbert space of a harmonic oscillator spanned 
by photon number states {\n) : n = 0,1,...}. Let N 
be the number operator acting as N\i) — i\i). For the 
Hadamard walk on the A^-circle we will employ a trun- 
cated Hilbert space Hn spanned by the N states with 
less than N photons. The random walk will take place 
on the phase states {|#o), • ■ ■ , \0n-i)} with 

N-l 

l^> = 7^E e ^ fe | fc >- ( 44 ) 

^ fc=0 

A shift on the phase states can be induced by the operator 
F± = exp(±^A r ), i.e. F±\6j) = \6 j±1 ). As in Sec. Eliot 
the two atomic states we employ be | f) and | }), and let 
S z be as in Eq. J3J. Then the conditional shift S is given 
by S = exp(-2iS z <g) ^-N) very similar to Eq. ©. To 
implement S we use a two-level model including ac-Stark 
shifts. The atomic levels are highly detuned from the 
cavity field, and the Hamiltonian for this effect is given 
by H = T%xN <£> 2S Z . The Hadamard transformation for 
the coin-flip on the atomic levels can be realized by a n/2 
pulse. The time r between successive pulses is set to be 
t = 2i:i/hxN. This implements the random walk. 

The standard initial condition would be to have the 
field in the state that projects to \9q) in the truncated 
space. This is unfeasible experimentally at this point. 
To make the random walk readily implcmentable with 
present state technology the initial state of the cavity is 
taken to be a coherent state \a) (a real and positive). It 
is possible to find such states \a) that have high overlap 
with \0q) and very small overlap with \0j) for j ^- 0, at 
least for N not too large |43|. 

Measuring the position of the walker corresponds to 
detecting a phase shift of an initial coherent cavity field. 
This can be achieved using homodyne detection (see 
for definitions). All ingredients employed in such an im- 
plementation of the random walk are achievable with cur- 
rent technology. 

Another upshot of this proposal is that controlled de- 
coherence can be introduced into the cavity. This allows 



to study the decrease in phase diffusion (transition from 
a linear spread a ~ T in the quantum case to a classical 
spread a ~ \JT) with noise (see also Sec. I Vlfl . Note how 
this decrease in dissipation contrasts sharply with the in- 
tuition provided by the fluctuation-dissipation theorem: 
that the introduction of noise usually results in increased 
dissipation! 

An implementation of random walks in optical lattices 
has been proposed by Diir et al. 0] • This proposal uses 
neutral atoms trapped in an optical lattice. The walk 
(on line or circle) takes place in position space given by 
the minima of an optical lattice. 

An optical lattice consists of periodic optical poten- 
tials formed by the standing waves resulting from two 
counter-propagating waves from two lasers. Here two 
identical one-dimensional optical lattices will be used, 
each of them trapping one of the internal states | t)> I I) 
of a neutral atom. The internal states can be hyper- 
fine states of the nuclear spin of the atom. The relative 
spacing of the potential minima of the optical lattices 
can be shifted by modulating the angle between the two 
counter-propagating waves, in particular the minima of 
the two lattices coincide when the angle is a multiple of 
7r/2. To implement the conditional shift we modulate 
the angle such that the lattice trapping the | t)-state 
moves with constant velocity v to the right, whereas the 
second lattice (trapping | 1)) moves with constant veloc- 
ity -ti to the left. This motion realizes the controlled 
shift. Laser pulses allow one to manipulate the internal 
state of the atom and to implement the Hadamard trans- 
form. For the random walk on the line the two lattices 
are initially tuned such that their minima overlap and the 
atom is placed into one of them, xo. The lattices are put 
into motion and the Hadamard pulse is applied at times 
t = nd/v, where d is the lattice spacing. A simple fluo- 
rescence measurement - together with several repetitions 
of the experiment - allows one to measure the resulting 
probability distribution. It is possible to implement a 
walk on the line with one or two absorbing boundaries 
by shining a laser at a certain position b which couples 
both internal atomic states to a rapidly decaying level. 
Modifications of the trapping geometry allow to imple- 
ment the walk on a circle |44| . 

This set-up is simple enough that many (hundreds) 
steps of the random walk could be implemented with 
present day technology. In particular since each lattice 
site is subject to the same coin-transformation, the pro- 
cedure does not require individual adressability of the 
different lattice sites, which is the big challenge when 
this type of implementation is used to build a quantum 
computer. 

Hopefully in the near future cither of this proposals 
or a new one will be implemented to begin experimental 
observations of quantum walks. Furthermore it would 
be interesting to see an experimental proposal for the 
continuous time random walk in some suitable physical 
system. 
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VI. DECOHERENCE 

The world of quantum mechanics shows all kinds of 
counter intuitive signs. Yet in the classical world of our 
daily experience we rarely see a cat that is in a super- 
position of dead and alive. The emergence of the classi- 
cal world from an underlying quantum micro-world is a 
topic of active debate, research and speculation. Usually 
it is assumed that interaction with a larger environment 
induces decoherence, which in turn makes the world clas- 
sical. We will describe how one could imagine that the 
classical random walk emerges form its quantum coun- 
terpart. 

The crucial difference between the quantum and the 
classical walk are the quantum coherences, which are 
present only in the first. So a transition from quantum to 
classical has to somehow kill these coherences. We have 
seen in Sec. IIIII that the classical walk can be obtained 
from the quantum walk if we measure the coin register 
at every iteration. Another way to obtain the classical 
walk is to measure the positions at every step. We can 
now imagine some continuous way to interpolate between 
non-measuring (the quantum case) and completely mea- 
suring (the classical case) and study the transition along 
this pathway. An intermediate step would be to measure 
the register (coin or position or both) only with some 
probability p mea s at each step, where p m eas = cor- 
responds to the purely quantum and p m eas = 1 to the 
purely classical case. Numerical studies of this transition 
from quantum to classical have been executed for the cir- 
cle, line and hypercube by Kendon and Tregenna and 
extensive analytical expressions for the walk on the line 
with decoherence of the coin-space have been obtained 
by Brun et al. g^lzl- 

They find that the classical behavior sets in very 
rapidly. As an indicator we can look at the variance of 
the walk. In the purely quantum case we have seen that 
a 2 ~ T 2 , whereas in the classical case a 2 ~ T. For the 
walk on the line with coin decoherence Brun et al. |49tl5C| 
compute the variance as a function of the transition pa- 
rameter Pmeas to find that even for very small values of 
Pmeas (very 'close' to quantum) the variance is already 
linear. Kendon's numerical results for both coin and po- 
sition decoherence confirm that for all p m eas > 1/T the 
behavior of the variance is essentially classical (linear in 
T). One could have expected that the variance would 
level off gradually with increasing p m eas, but this does 
not seem to be the case. 

There is another interesting option to transition the 
coined walk between quantum and classical taken by 
Brun et al. 0, |mJ . Imagine the quantum discrete time 
random walk with the modification that at each iteration 



we use a fresh coin Hilbert space to flip our coin. This 
does not allow the quantum coherences between the coin 
degrees of freedom to ever interfere and hence the walk 
behaves like the classical walk. A way to go from quan- 
tum to classical would be to use a new coin space only 
once in a while; a walk alternating between two coin- 
spaces at every step can be thought of close to the quan- 
tum walk, a walk with three coin-spaces is already less 
quantum and so on to the walk with an amount of coin 
spaces close to T, which is close to classical. It turns 
out that this way to transition from quantum to classical 
induces a very different behavior of the variance a 2 : it 
stays quadratic in T for a whole range of coin-spaces to 
become linear only in the limit of one new coin per time 
step. So the quantum character in this transition is very 
robust! 

Kendon and Tregenna's |48j | numerical studies also 
show that decoherence can improve both convergence to 
uniform and reliability of hitting. For small p <C 1/T 
(where the walk is already a 'little bit classical') under 
certain conditions the convergence to uniform on the cir- 
cle (or line-segment) can be faster than in the purely 
quantum case. 

There is certainly a lot more to understand about deco- 
herence in quantum random walks than we do presently 
- these initial studies seem to unveil just the tip of an 
enormous iceberg. 

Conclusion: Starting from a rather simple physical ef- 
fect we have tried to build up the beautiful framework 
of quantum random walks. We have presented several 
rather striking differences in the behavior of these walks 
as compared to their classical counterparts concerning 
their distribution, their mixing and hitting properties. 
With the advent of quantum computation it is our hope 
that quantum random walks will be utilized to provide 
new and fast algorithms that can be run on a quan- 
tum computer. We have seen that the first steps have 
been made in this direction, also from the implementa- 
tion point of view. This area of research is still dotted 
with open questions and problems and we hope the reader 
has obtained a good general intuition and overview as to 
how this exciting field has evolved so far and where it 
might go. 
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